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Abstract. 

Instanton was found for Kraichnan model with 'frozen' velocity field. Large order 
asymptotic of quantum-field perturbation expansion for renormalization constant Z v 
was investigated. It was shown that this expansion is convergent one. The radius of 
convergence was calculated. 



1. Introduction 

Large-order asymptotic analysis of quantum-field perturbation expansions is an 
actual problem of modern statistical physics. Direct perturbation calculations are 
cumbersome and difficult. Knowledge of large-order asymptotic and choice of the 
corresponding resummation procedure allow to get more or less good estimation on 
the base of a few first terms of perturbation series. Application of the resummation 
procedure without knowledge about large-order asymptotic behaviour can lead to 
inaccurate results. 

Article [1] can be considered as a first attempt to solve this problem. The 
asymptotes were estimated here merely by number of graphs in the perturbation 
expansion order. The correct large-order asymptotic investigation proposed by Lipatov 
[2] is based on the saddle-point calculation of path integrals. It was used for all main 
quantum-field theory models and static models of critical behaviour, see [3]. Large- 
order analysis also was consistently constructed for the dynamic models with equilibrium 
static limit [4, 5] using standard Martin-Siggia-Rose (MSR) [6] variables. For all models 
mentioned above the instanton was found and divergent character of perturbation series 
was proved. Then there is common opinion that the instanton existence always leads 
to divergent series. 
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This paper accounts the Kraichnan model with 'frozen' velocity field. It describes 
the turbulent diffusion in stationary random field [7, 8, 9] and well known problem 
of random walks in random media [10] [15]. The renormalization group (RG) 
approach was used for investigation of the scaling behavior in this model [8, 9]. 
The objects of interest are renormalization constants. In this paper large order 
asymptotic of perturbation expansion of viscosity renormalization constant in this model 
is investigated on the base of instanton found and it is shown that the coefficients of 
perturbation series grow essentially slower than N\. The series of perturbation theory 
turn out to be convergent. The radius of convergence is calculated. 

Then our result contradicts the common opinion about connection between 
instanton existence and perturbation expansion divergence. The more or the less general 
explanations of this phenomenon is given in Sec. 4. Let us note that the examples 
of instanton analysis of convergent series are known. One of these is the convergent 
perturbation series introduced in [16]. The instanton analysis of the convergence of this 
expansion was fulfilled in [17]. Another example is standard Kraichnan model [18, 19]. 

The large order asymptotic analysis is more difficult in dynamic models then in 
static ones. Usually there is no instanton in natural class of MSR variables here. 
For example the absence of instanton in Kraichnan model of turbulent diffusion was 
proven in [20]. Fortunately, the Lagrange variables can be used in this model [21, 22]. 
In these variables instanton was found and the large-order asymptotic behaviour was 
investigated [22]. The perturbation series in this model appeared to have a finite 
radius of convergence. Note that a specific feature of standard Kraichnan model is 
a proportionality of the velocity field correlator to 5(i)-function. Then a lot of graphs 
of perturbation theory are absent here [23] that could explain the convergence of series 
discussed. 

Because of the difficulties of instanton analysis the large-order asymptotic form was 
estimated in some papers merely by the number of graphs at large order of perturbation 
expansion in spirit of [1]. This approach produces accurate results for models with 
scalar and vector fields without derivations in interaction. But we will show that such 
estimation [7] may cause a mistake. Kraichnan model with frozen velocity field discussed 
in this article is a good example of this fact. In contrast with standard Kraichnan 
model, here the velocity field correlator doesn't depend on time, so the number of 
perturbation diagrams demonstrates a factorial behaviour N\ while the perturbation 
series is convergent. Simplified example of model considered with the constant velocity 
was considered in [24], where the number of perturbation diagrams demonstrates a 
factorial behaviour too, the instanton was found and the perturbation series convergence 
was proved by exact solution of the model. 

This paper is organized as follows. The Kraichnan model with a "frozen" velocity 
field is described briefly in Sec. 2. Also the response function to be studied using 
MSR-formalism is introduced. The composite operator in Lagrange variables which is 
used for calculation of renormalization constant Z u is introduced in Sec. 3. Instanton 
approach for this composite operator is presented in Sec. 4. A particular solution 
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of stationary equations used in the following analyses is investigated in Sec. 5. The 
renormalization of the Green function with the composite operator is described in Sec. 
6. Large-order asymptotic for expansion of In Z u is calculated in Sec. 7 using the replica 
trick. Cumbersome stationarity equations for some numerical parameters are discussed 
in Appendix. 

2. Kraichnan model with a 'frozen' velocity field 

Kraichnan model describes the turbulent advection of passive scalar admixture in 
rf-dimensional fluid. It is based on a stochastic equation 

[d t + gViVii*, t) - z/A]y9(x, t) = £(x, t). (1) 

Here x6l d and t are space and time variables, </?(x, t) is a passive scalar field, V(x, t) 
is a random vector velocity field, £ (x, t) is a random force, v is a viscosity, g is a coupling 
constant. Laplacian and gradient refer to the space variable x; here and henceforth all 
derivatives in squared brackets act on the field tp as well. For shortness we introduce 
d t = d/dt. The convolution with respect to repeating subscripts is implied here and 
henceforth. 

The random values £ and V are supposed to be distributed by Gauss low. It is 
known that the results obtained in RG analysis are independent from [8, 9] then 
one can state the correlator has an arbitrary form. In contrast with the standard 
Kraichnan model in the model discussed the velocity field correlator doesn't depend on 
time 

(V i (x,*)V j (x',0) = A,(x-x') 

(compare with (Vj(x, i)Vj(x', t')) = Aj(x — x.')5(t — t') in the standard Kraichnan 
model), in other words one can consider V-field as a time independent. The velocity 
field correlator in the momentum representation has the power-like form [8, 9] 

Dg(q) = / M 3 (z) exp(.qz) = A T (<5« - ^) ^ + A,^^, (2) 
that was used for RG-analyses of the model. Then in the coordinate representation 

D tJ (z) = a x -| + a 2 ^, P = d/2-a = l- e/2, (3) 
where the parameters a±, a 2 are known: 

^- 22Q+ Jff (a + 1) (M2^-1) + A L ), (4) 

_ (\ x ^ r(^+i) ( . 

^-(At-A l ) 22and/2r{a + iy (5) 

At and are transverse and longitudinal coupling constants. The parameter A^ 
corresponds to compressibility of fluid. 

The model considered is important for the description of diffusion in random fluids 
[8, 9]. Moreover it relates to the developed turbulence problem [25]. In fact the model 
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(1, 2) has two coupling constants Xt, Xl- Nevertheless for large order asymptote of 
e expansion investigation these can be reduced to the only one coupling constant g. 
Indeed in a fixed point all coupling constants are proportional to a small parameter 
of regular expansion e = 2 + 2a — d, [8, 9] Xt ~ ~ g 2 . The similar situation was 
observed in dynamic models C-H (A, B, C,... are a common nomination for particular 
dynamic models introduced in [26]) with equilibrium tp A static limit [27, 5]. 

The infrared behaviour of the model was investigated by means of RG method. 
The renormalization yields the substitution v — > vZ v , g — > gZ g , the renormalization 
constant Z v is calculated by means of perturbation theory, then it has a form of series 
in coupling constant. The properties of this series are the main subject of this paper. 

Usually MSR-formalism [6, 27] is used to transform stochastic models to quantum- 
field ones. The base equation (1) is represented in a form of path integration in 
auxiliary field tp'; path integrations in tp and V fields are introduced to study averaged 
characteristics of fluid. Then the expression for the response function in an arbitrary 
velocity field has a form 

r (i,2) / WW ?(xi, *iMx 2 , t 2 ) exp(S msr ) 

' f WW eMS\?=o) ' { ) 

with the renormalized action [8, 9] 

S msr = V?/(X ' t)jP 2 ^ /(x/ ' t/) + p'(*,t)[d t + gZgViVifr) - uZ v A]<p(x,tU7) 

usual standard agreements for dynamic models [27] and all integrations needed are 
implied henceforth. The normalization factor in Exp. (6) corresponds to a free model 
(at g = 0). 

After the integration in V field one obtains MSR representation for the renormalized 
response function 

f V VG (1 ' 2) e' v * (x)D ^ 1 (x ~ x,) v ^ x ')/ 2 



3. Lagrange variables 

As it was stated in papers [22, 20] there is no instanton for Kraichnan model in the 
framework of MSR formalism. Similar arguments are correct for the model with the 
frozen velocity field. But Lagrange variables [21, 22] can be used for instanton analysis 
of the model considered as in standard Kraichnan model. 

Seems these variables have an origin in the quantum-fields methods application 
in random walks and macromolecules problems (see [28]) in analogy with Hamiltonian 
form of the standard Feynman-Kac path integral. Lagrange variables can be introduced 
successfully in the dynamic models with linear in main field tp stochastic equation (1) 
only. Then the dynamic equation for the response function (6) can be considered as 
Schrodinger equation [21] or as Fokker-Plank equation [22], and response function in an 
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arbitrary velocity field can be represented in the form 

r (i, 2) = e(t 2 -tx) fDcTtcf eMS L9r ) . , 

v " (47n/) d / 2 (t 2 -ti) d / 2 JDcDc' ex V (S\ L g g r oY 
ti 

S L9r = J dr(^ - uZ u c'\r) + ic\r)d r c{r) + g Z g c' (r)V (c(r)^j . 
The boundary conditions 

c(ti)=xi, c(t 2 )=x 2 

are implied for the path integration in c field in the numerator. The integration in c' 
fields is supposed to have free boundary conditions. The normalization path integral 
corresponds to a free model with g = and may be calculated at zero boundary 
conditions for c, c' fields. Vector fields c(r), c'(r) play a role of coordinates and momenta 
of fluid particles and depend on time only. 

Let us note that this representation produces one Green function of the model 
only, namely the response Green function in an arbitrary velocity field. Then the MSR 
action (7) can be obtained from (9) by no change of variables and the statement about 
an instanton absence in MSR variables is not correct for Lagrange ones. 

Lagrange variables can be used to investigate renormalization constant Z u in theory 
(7). Indeed let's differentiate Exp. (8) with respect to v in order to extract Z v constant. 
Then one has two point Green function with the (p'vZ u A(p composite operator insertion. 
Using the response function (8) this Green function 

Z u < (p(xi,ti)ip'(x 2 ,t 2 ) J dx dt Lp'(xo,t )A(p(x ,t ) >= G 

can be rewritten as 

G = Z v J dx dt (10) 

/ !DVGy' 2 \xi, ti, xp, t ) AG™ (x , t , x 2i t^MW 

/ DVe" Vl(x)D ^ 1(x ~ x ' )Vj{x ' )/2 
The renormalization constant Z v has poles in e at e — > that have to cancel UV 
divergences of the model. As we consider a renormalized response function (8) Exp. (10) 
must be finite at e — > 0. Then 

res \nZ u = — res In / dx. dt G (11) 

e^O e^O J 

and the Green function G contains all information needed about the poles of 
renormalization constant Z v . 

The diagrams for G include the internal loop part and the external full propagator 
part without divergences. The last does not contribute into Exp. (11) then one will 
discuss now the amputated diagrams for G. 
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The Green function G can be easily written in Lagrange variables. The following 
Gaussian path integration in field V and Fourier transforms in x 2 — x x and t 2 — t\ 
variables yield the action to be studied 

S = -iq(x 2 - xi) - vZ v (c? + c'i) + ic^dd + ic' 2 dc 2 + Z U S U , (12) 

where q is a momentum. Frequency u = 0, that is sufficient for the renormalization 
constant investigation, because this constant is frequency independent. A nonlinear part 
of the action is collected in the term 

Su = -f (^(tOA.MtO - c 1 (r[)y ij (r[) + 

+4(r 2 )Aj(c2(r 2 ) - c 2 (r 2 ))c 2j (r 2 ) + 2c' H Ai(ci - c 2 )c 2i ^ , u = g 2 . 

Here and henceforth one implies that the fields c u c[ (I — 1,2) with argument omitted 
depend on T\ . All necessary integrations in T\ and the ranges of integration 

h < n, t[ <t < r 2 , r 2 < t 2 

are assumed. Finally one gets for Fourier transformed G function the following 
expression in Lagrange variables: 

G= I d(x 2 -xi) / dih-h) 1 (13) 

Cl(tl)=Xl, C2(*2)=X2, 

J 'Dc 1 Dc 2 'Dc' 1 Dc' 2 We s , 

Cl(to)=C2(to)=XO 

where fore-exponential factor 

/ dTidT 2 (ic[ + (ti - r 2 )9 r i) / dT X dT 2 {ic! 2 - (t 2 - r 2 )9 r i) 



W 



(h - t )(t 2 - 1 ) 

j rfrirfr 2 (ti - 7i)(t 2 - r 2 )9 r 2 



(14) 



^ = uZ u c' u — Dij(ci - c 2 )c 2j 



(ti-t )(t 2 -t ) 

'del 

is produced by Laplace operator in (10). Let us note that the non-trivial boundary 
conditions for Ci, c 2 fields 

Ci(*i) = xi, c 2 (t 2 ) = x 2 , d(to) = c 2 (t ) = x (15) 

mean the path integration with fixed boundary conditions. The integration in c[ fields 
is supposed to have free boundary conditions. The normalization factor K corresponds 
to a free model with u — and may be calculated at zero boundary conditions for cf. 

Cl(tl)=C 2 (*2)=0 

K= J DciDcaDc'iDc'a exp(5| 9=0 ). (16) 

ci(to)=c 2 (to)=0 
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For shortness we introduce the parameters 

T = t 2 - tl , Tx^fo-ti, T 2 = t 2 -t , (17) 
x = x 2 - x x , x (1) =x -xi, x (2) =x 2 -x , (18) 
the values x, xo, x±, x 2 , x^\ x^ are the modules of corresponding vectors. 

4. Instanton analysis 

In order to extract the iV-th term of perturbation series the Cauchy formula is 
traditionally used [2]: 

° m = hflM dU - GW = EG'V, (19) 

J N=0 

the integration is produced along a closed contour containing zero in a complex plane. 
Let us extract the large N parameter from the action S (12) by the following scaling 

{c i)C ;} -> {N^c^N 1 / 2 ^}, x^v^x, u^N^u, (20) 

the scaling of x variable is necessary due to the connection between c and x based on 
(15). The same scaling in normalization factor K cancels the determinant corresponding 
to this change of variables. The momentum q is scaled q — > iV 1//2 q also. It is possible 
because the renormalization constant investigated is momentum independent. 

The integrals in c', c and u at large N can be calculated by the saddle-point 
approach. The main contribution at N — > oo is given by the integration near the 
instanton that is a special realization of variables c st , c' st , u st . The action has an 
extremum at the instanton. 

Let us mark that the scaling procedure is an essential in the determination of the 
connection between instanton existence and divergence or convergence of perturbation 
theory. Usually the quantum- field theory action has a form S = Sq + Si n t, where Sq is 
a free part of the action, So — where $ denotes a field or a set of fields of the 

theory considered, K is a linear operator. Si n t is an interaction of the form Si n t = A$ fe 
(k > 2) in a local theory. Presence of derivations in the interaction does not affect 
analysis presented below. A is a coupling constant and the expansion parameter. To 
extract large iV parameter from the action the scaling $ — > \/N§, X — > X/N^ 2 ^ 2 is 
necessary. Then the main exponential contribution in expression similar to (19) due to 
instanton is proportional to 

where S st is an action in the stationary point. This contribution demonstrates the 
divergence of the perturbation expansion in A due to k > 2. Then the divergence of the 
series is connected with the scaling of the coupling constant. In [3] another instanton 
analysis scheme without Cauchy formula using was proposed. But the results of analysis 
presented above will be the same in both schemes. 
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Note that the model considered (12) in Lagrangian variables differs from the general 
case. The coupling constant is scaled here by positive power of N due to the nonlocal 
character of the interaction, then instanton analysis can lead to convergent series as in 
standard Kraichnan model [22] or simplified example of Kraichnan model with constant 
velocity field [24] , where the instanton was found and the perturbation series convergence 
was proved by exact solution of the model. 

Let us describe instanton calculation in the model (12) in more details. It's quite 
reasonable to simplify problem by taking into account the symmetry of the model that is 
initially violated by x vector only. Then we suppose that the fields q, cj are parallel to x 
and the only modules q, c\ must be found. Let us mark that the stationarity equations 
are non-linear differential ones. The existence and the uniqueness of the solution are not 
proved in general case. We propose to use the solution with the same symmetry as Green 
function investigated. Note that the spherical symmetry of the instanton in Lipatov 
work [2] based on the same ideas. The possibility of other solutions existence with 
other contributions to large order asymptotes is an open question. The same situation 
is observed in every case of instanton analysis. The supposition about the symmetry 
of solution used here was proven in [24] for simplified Kraichnan model with known 
exact solution. For standard Kraichnan model considered in [22] this supposition yields 
the result coinciding with the exact known anomalous dimensions of a set of composite 
operators. 

All stationarity equations are supposed then to be projected on the direction of x 
vector. It also simplifies the tensor structure of D correlator that has more compact 
form now: 

£(x) = ^p D = ai + a 2 . (21) 

Except the regular in e terms the action S (12) contains poles in e due to Z v , Z g 
constants. It was shown in [29] that while a renormalization constant is investigated 
the corresponding singularities must be extracted before any instanton calculations and 
later they contribute only in a fore-exponential factor of the saddle-point method. This 
extraction of singularities is necessitated by existence of two large parameters, namely 
1/e connected with regularization and the saddle-point method parameter N. Due 
to the renormalization approach the value Ne must be considered as a small one in 
the framework of instanton analyses [29]. Then the exponential term in (13) must be 
presented in a form 

OO j 

exp(S) = exp(S reg + S sing ) = exp(S reg ) ^ — ( S sing) p , (22) 



pi 



<7 

'-'reg 



Z u =l 
Z v =\ 



S smg = -v{Z v - l)(cf + cf) + (Z u - 1)S V 



and only the term S reg in the exponent must be variated. In standard Kraichnan 
model this approach was proven in [22] by the comparison of the radius of convergence 
calculated with the exact known results. Mention should be made that the l.h.s and the 
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r.h.s of the identity (22) are essentially different under the integral § du/u N+1 and they 
yield different results of the saddle-point method due to the competing of parameters 
N and 1/e. 

Let's remind the renormalization constants Z u , Z u in minimal subtraction (MS) 
scheme have a form 1+ poles in e terms, so that [Z — 1) contain the pure singularities at 
e — > only. Exp. (22) shows that the divergences in e make a sense in the framework of 
perturbation theory and diagrammatic expansion only. Then the path integration must 
be interpreted as a sum of perturbation terms, the renormalization is supposed to be a 
cancellation of divergences. 

Then let us set Z v — 1, Z g — 1 in the expression for the action (12) in order to 
write the regular instanton equations. The variations of action S involve the integral 
operators of the form 

[Afcc' fc ](C) = j dr k D( Cl (C) - c fc K, l,k = 1,2. 

For example the variation in c x yields 

=► uc'^Odt ([D u c5i](0 + [^i2c' 2 ](C)) = -id^d^m 



MO 

as well the variation in c[ 

J^L = =► -2z/c' 1 (C)+^ f ci- M ([ J D 11 c' 1 ](C) + [ J D 12 c / 2 ](C)) = 0.(24) 

The contribution of integral operators in eq. (23) can be excluded with the help of 
eq. (24). For this purpose we should differentiate (23) in ( 

d( 6^) = ° ^ -2^c c 'i + " u $t ([^n c i] + [^124]) = 0, 

express [DncJ] + [Di 2 c' 2 ] and substitute it in the eq. (23). It yields the 2nd order 
differential instanton equation 

- 2vd c c' 1 + t 9 -4j^ = (25) 

that can be solved with respect to d^c(Q. The same calculation can be done for c 2 field. 
The solution of equation (25) has a form 

0^(0 = ^-^(0, 1 = 1,2 

and contains the arbitrary parameter F. Each numeric value of F constant corresponds 
to a particular solution with its own boundary conditions. 

The following calculations at arbitrary F can not be performed analytically, the 
simplest case at F = allows to reduce the instanton to quadrature only. Nevertheless 
the corresponding particular solution can be used to determine the asymptotic behaviour 
of the renormalization constant investigated. 
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5. Particular solution 

Let us consider the particular solution with F = and determine its boundary 
conditions. The solution 

cJ(C) = ^, 1 = 1,2 (26) 

can be substituted into the variation equations with respect to c\. Using the identity 
dndci = dQ the result can be written in a form 

X2 



-d& = l J D( Cl (C)-z)dz, 1 = 1,2. 



XI 



The last differential equation can be easily integrated and this leads to the solution for 
the q(C) fields in quadrature: 

ci(C) c 2 (C) 

dc _ u(( - ti) f dc _ u(( - 1 ) ^ 



i:;D(c-z)dz v ' J %D{c-z)d 

XI x 

Note that an analytic regularization is assumed in (27) that eliminates the singularity 
in z = c point. 

After the substitution of explicit form D(c) (21) Exp. (27) produces the boundary 
condition interested 

xW/x x^/x 

— f f( v )dv = - UD ° r = — [ f(v)dv, (28) 
Ti J K J x 2 - £ {l-e)v T 2 J K J 



where the function f{v) introduced is 

v'-Hl-v) 1 -" 

f(v) 



v i-e + (i - v y- £ ' 

One sees that the case F = is simple enough to give a quadrature representation 

(27) for ci, c 2 fields with the boundary conditions (15). 

Initially the problem consists in calculation of instanton for an arbitrary boundary 
condition x(T). Exp. (28) solves the problem in the specific case with a special value of 

x^ 2 \ and x. Then one has instanton for functional integral in c, c' fields. But the 
object investigated (19, 12) includes integrations in x, T, x , t and u also. Our main 
idea here is to explore the independence of the renormalization constant investigated on 
the momentum q. Let us include the integrals in variables x, x , to into the saddle-point 
method. Let's choose the momentum q so that the solution of a stationarity equation for 
x variable be exactly equal to the result obtained with the help of boundary condition 

(28) corresponding to the case F — 0. This choice gives us a chance to explore the 
particular solution constructed analytically and to solve the problem without numerical 
calculations. Moreover constant F is not a free parameter in this approach, then one 
has no problem of zero modes connected with F arbitrariness. 
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The stationarity equations for x, x , to variables are too cumbersome to be written 
down here. Nevertheless, as shown in Appendix of the article, these can be simplified 
significantly using particular solution (26), its properties, and integration by parts. Then 
the equations solution calculated corresponds to the case = x^ 2 \ 7\ = T 2 and the 
action has the following form in stationary point (26) 

g rt = -tqx- 2 ^f° ^ . (29) 
v l e[\ — e) 

Now the stationary equations for x in F = case: 

2 X 

J D(z)dz. (30) 



SS n iqv 2 
— = =s> 
ox u 





Besides we have the boundary conditions (28) imposed by our choice of particular 
solution F — 0. By solving equation (30) one obtains the proper value for q in case 
F = 0: 

a - an - iD ° u x 2 ~ £ - x 2 ~ £ - UD ° T/2 (31) 

9-90- (1 _ £ wi- £i/2 ' x ~ x st ~ 1/2 • W 

(l-e)u J f(v)dv 
o 

Combining (29) and (31) finally, one obtains the action S (12) at the stationarity solution 
and Z u — 1, Z u — 1: 



1st 



v 2 e 



(32) 



6. Simple poles in e 



Due to (11) simple poles in e of G contain all the necessary information about the 
poles of renormalization constant Z v and the corresponding critical indices. 

Since Dq(e) = a± + a 2 (21) and (4,5) Dq(e) function can be presented in a form 

D (e) = A + eB(e), B(e) = B + B x e + 0(e 2 ). (33) 

It now follows that the action (32) has a form 

_ uA uA{x £ st - 1) uB{e)x £ st 

&st — o 2 2 ' I 

The first term here is singular in e — > 0, then as well as S S i ng term discussed in 
Section 4 it must be presented in fore-exponent form (22). The second term seems to 
be finite at small e. Nevertheless its logarithmic behaviour in x st results in singular 
in e contribution to the large order asymptote. We will discuss this at the end of this 
Section. As a result the only regular term of action S st (34) is the third one. Then Exp. 
(22) must be corrected by changing S reg -> S reg , S sing -> S sing , 



Sainn $ sir. 



uA uA(x e st - 1) 



'smg '-'sing V 2 £ 
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Combining Exp. (34) with the result for x st (31) we get 
S reg = (uT^r /(2 - £) P(e), P(e) = ( 



B(e) ( D, ^ e/(2 - £) 



e)v2jl /2 f{v)d 



Thus the Green function studied is of the form 



G [N] = N -N(i-e) j _du_ j dTZ ( TjUj£)eMNSreg (e,T))x 



oo 1 

x EriW( 1+ °( r1 ))- ( 35 ) 

p=0 

The Z factor stays for a Gaussian fluctuations contribution and fore-exponential 
factor W (14), the corresponding integration is normalized by K factor. 0(N~ 1 ) term 
shows the accuracy of calculation at iV — > oo. 

The amputated Green function with composite operator considered must be 
dimensionless. The factor T~ l restoring this zero dimension has to be produced by 
Z and we will extract T _1 from Z in order to stress this fact. 

The integration in T then diverges as a logarithm and produces singularities. It can 
not be treated by the saddle-point approach. The analogous situation exists in a well- 
developed instanton analyses for static ip 4 model where the role of divergent integration 
parameter T plays the scale parameter in the coordinate space [2, 29]. 

Let us change variables 

u = uT £ '\ 

the new variable u is dimensionless. As a result the expression 



x 



x | 



x LriW p ( 1 + °( r1 )) ( 36 ) 

p=0 

can be integrated over T in UV region (small T). This yields a simple pole 2/ (Ne). The 
convergence of integral at large T is provided by IR regularization assumed. Note the 
factor Z does not contribute to the stationary equations as it does not depend on N. In 
the same way this factor did not affect the simple pole in Ne at least at principal order 
in 1/N. 

The integration in u is investigated by the saddle-point approach. The stationary 
equation with respect to u is 



d 



2 



P(e)u 2 -c — In if 



0, then u 



2 
2-e 



1 -e/2 



du • st P{e) 

So we obtain the leading order in N 

2 {e/2 _ lU P(e)e x(iv + i)(i- £ /2)^ i f _ 

p=0 
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The factor C(e)N p appears due to Z(u st , e) contribution and the fluctuation integration 
in u; p is a constant. 

Let us discuss the residue in e calculation. Simple poles in e can appear if higher 
poles of S S i n g in the fore-exponent are multiplied by regular in Ne contribution of the 
exponent term exp(S reg (e)). Therefore in the MS scheme chosen all p 7^ terms of the 
sum contribute to the result. Fortunately a finite renormalization could help us to map 
out all this terms. Indeed let's scale the expansion parameter 

uD (e) — > x(e)«D (e). 

Such a renormalization is equivalent to a scaling of the velocity field correlator and 
doesn't affect the scaling dimensions. Let's choose k(e) so that the regular part of 
action loses its dependence on Ne, namely 



= K . 

e=0 



f f{e)e \ / f\e)e \u- £ 

Vl-e/2/ ~* \\-ej2) 

The corresponding renormalization can be easily written as a perturbation series 

x{e) = l+(—M-^f\-^-V + 0(e 3 ), 
v ; \2A L 6Au\ A J K h 

the parameters A, B , B x are introduced in (33). After the scaling all poles in e contained 
in the sum with respect to p don't contribute to the simple pole. As a result the residue 
discussed demonstrates the asymptotic behaviour 

res G [N] = Const N Const K^/N\, N -> oo 

that corresponds to a finite radius of convergence for the perturbation series of G 
function. 

Now let's show that the second term in (34) doesn't change our answer. Indeed 
its behaviour in x is logarithmic. Then the presence of additional factor \nx ~ ln(T) 
in (36) is equivalent to the additional operation d/d(Ne) of T Ne ^ 2 factor in (36). This 
operation can't produce simple pole in e as our expression doesn't contain logarithmic 
in e contributions. It was shown in details for the similar problem in [29] . 

The following step is to calculate \nG(u), extract the residue in e = of the N- 
th order of perturbation theory and explore formula (11). It's useful to transform the 
logarithm of the composite operator G with the help of replica trick. 

7. The replica trick 

The simplest way to present a logarithm of arbitrary integral expression in an 
integral form is to use the formula 

In f dTf(T) = lim |- JJ / dT a f(T a ). (37) 
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As a result the variable T becomes an r-dimensional vector in a replica space [27]. 

Performing such a procedure with respect to G one gets the expression similar to 
(35) where the integrals must be rewritten as follows 



r— 1 

f^Ti J(ndT a )eMNS reg )Z({T a y-J ), 



(38) 

Sre 9 = XW 2 ) 2/(2 ~ £)P ^ 
o=0 

the factor Z depends on the replica variables {T a }. Nevertheless one knows that the 
expression investigated with the help of replica trick is also dimensionless due to Z 
contribution. Besides the expression must be proportional to r at small r. Then the 
operation limd/dr in (37) yields a non-trivial finite result. 

r^O 

It is easy to see that the saddle-point approach can not be applied to all integrations 
in (38). At least one of integration considered has a non-saddle-point structure. We 
observed the same situation in the previous section in case of integral in T. 

Let us exclude the integration in To from the saddle-point approach consideration. 
The set of other variables {Tq,}^~\ has a zero solution at the saddle point. So the 
stationarity solution for u s t depends on To only. As this non-saddle-point mode in the 
replica space is chosen in an arbitrary manner, in fact we have constructed r identically 
instanton solutions in the replica space. This results in the factor r which allows us to 
produce the operation \imd/dr correctly. 

r— >0 

The integration in T must be treated in the same way as the integration in T in 
the previous section and yields the same result 

res In Z u = ConstN Const K* /N\, N -> oo. (39) 

The calculation of Const in this formula could not be produced without an explicit 
calculation of the fluctuation integral. But this is difficult and cumbersome problem 
that is not to be solved here. Thus the asymptotic form (39) demonstrates that the 
series investigated has a finite radius of convergence. 

In fact this section can be resumed as follows. We have shown the saddle-point 
method yields the appropriate result namely the circle of convergence for function In G 
is determined by the same singularity as function G. The properties of this singularity 
were calculated by saddle-point method. Other singularities for InG could exist in 
principle but these have non-saddle-point structure. 



8. Conclusions 



We have constructed the family of instantons in Kraichnan model with "frozen" 
velocity field and found one of them in explicit form. Considering the asymptotic 
behaviour of the renormalization constant at large order of perturbation expansion we 
have demonstrated that the corresponding perturbation series has finite convergence 
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radius. Furthermore, we have disproved the common statement that the behaviour of 
the series may be defined by quantity of diagrams at large order of perturbation. 

Our results can be used for an improvement of resummation procedures constructed 
for the Kraichnan model. 

The article was supported by RFBR (Grant No.08-02-00125a) 2008-2011. 



9. Appendix 



Let us consider the stationarity equations for x, x , to variables. The linear 
translation of fields 



Cl (n) = x.+xW^-^+c^), C2 ( r2 ) = xo+x( 2 >^^ 



fc 2 (r 2 )(40) 



Ti iV A " ^ ^ " ' T 2 
shows explicitly the dependence of action S on Xi, x 2 . The boundary conditions for 
new fields q are assumed to be zero. The differentiation of action S in x ,to produces 
cumbersome terms due to the interaction part S u of the action. These terms are of the 
form 



u 



71 



(41) 



where D' is a derivative of correlator D on its argument (21). 

Due to D correlator is an even function, it's enough to calculate in Exp. (41) only 
the term corresponding to the T\ contribution. Substitution of particular solution c[ 
(26) yields 

dD{ Cl {r[) - ci(n)) 



U f A > dc ^ T 'l) f A 

v*Tl J ~dr[~ I 1 



dn 



The inner integral can be calculated by parts. Then the integral term is proportional 
to J rfr{[_D 11 c' 1 ](r{) and surface terms can be calculated trivially. The result turns out to 
be calculated with the help of the following identities based on the stationarity equations 



d( 



5S 







d(u[D ll c' l \(0+u[D 12 (c' 1 + c' 2 -c' l )}(0 



I 



dC(-2i/cJ(C) + idccO 



I = 1,2. 



The last equality is written using the eqns. (26). 

In this way the stationarity equation discussed has the form simplified by F = 
condition: 

:r( 2 ) 

r (2) r (i)l r 



5t 



c (i) 



SS 
fa 



iqv 
u 



X 

I 



D{z)dz. 
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Besides we have the boundary conditions (28) imposed by our choice of particular 

solution F — 0. Let's note that the non-trivial instanton equation with respect to 

xo variable turn out to be an identity at F — 0. Indeed the substitution of (26) into the 

action S (12) yields S(x) as the function in the unique x = x 2 — Xi space variable. But 

the boundary conditions (28) give an oportunity to determine stationary value of x . 
[1] A.G.Basuyev, A.N.Vassiliev Sov. Phys. TMF, 18, p. 181 (1974). 
[2] Lipatov L. N., Sov. Phys. JETP, 45, p. 216 (1977). 

[3] Zinn- Justin J., Quantum Field Theory and Critical Phenomena. Oxford Univ. Press, Oxford, 
(1989). 

[4] J. Honkonen, M.Komarova, M.Nalimov, Nuclear Phys. B, 707 (FS) p. 493, (2005). 

[5] J. Honkonen, M.Komarova, M.Nalimov, Nuclear Phys. B, 714 (FS), p. 292, (2005). 

[6] P. C. Martin, E. D. Siggia, and H. A. Rose, Phys. Rev. A, 8, p. 423 (1973). 

[7] Steven A.Orszag and Victor Yakhot, Journal of Scientific Computing., 14, No. 2, p. 147 (1999). 

[8] Juha Honkonen and Esa Karjalainen, J. Phys. A:Math.Gen., 21, p. 4217, (1988). 

[9] Juha Honkonen and Esa Karjalainen, Physics Letters A., 129, p. 333 (1988). 
[10] Marinari E, Parisi G, Ruelle D and Windey P, Phys. Rev. Lett., 50, p. 1223 (1983). 
[11] Marinari E, Parisi G, Ruelle D and Windey P, Commun. Math. Phys., 89, p. 1 (1983). 
[12] Fisher D S, Phys. Rev. A, 30, p. 960 (1984). 

[13] Fisher D S Friedan D, Qiu Z, Shenker S J and Shenker S H, Phys. Rev. A, 31, p. 3841 (1985). 

[14] Kravtsov V E, Lerner I V, and Yudson VI, J. Phys. A:Math. Gen., 18, p. L703 (1985). 

[15] Kravtsov V E, Lerner I V, and Yudson VI, Zh. Eksp. Tcor. Fiz., 91, p. 569 (1986) [Sov. Phys. - 

JETP 64, p.336 (1986)]. 
[16] Ushvendze A. G. Yad. Fiz. 38, p. 798.(1983). 

[17] Honkonen J., Komarova M., Nalimov M. Acta Physica Slovaca, 52, p. 303, (2002). 

[18] A.M. Obukhov, Izv. Akad. Nauk SSSR, Ser. Geogr. Geofiz., 13, p. 58 (1949). 

[19] R.H Kraichnan, Phys. Fluids, 11, p. 945 (1968). 

[20] Balkovsky E., Lebedev V., Phys. Rev. B, 58, p. 5776 (1998). 

[21] M. Chertkov, Phys.Rcv. E, 55, p. 2722, (1997). 

[22] A.Yu. Andreanov, M.Komarova, M.Nalimov, J. Phys. A, 39, p. 7801, (2006). 
[23] Loran Ts. Adzhemyan, Nikolaj V. Antonov, Phys. Rev. E, 58, p. 7381 (1998). 
[24] M.V. Komarova, I.S. Kremnev, M.Yu. Nalimov, Theoretical and Mathematical Physics, 158, p. 167 
(2009). 

[25] J. P. Bouchaud, A Georges, Phys. Rep., 195, p. 127 (1990). 

[26] Hohenberg P.C., Halpenn B.I., Rev. Mod. Phys., 49, p.435 (1977). 

[27] Vasiljev A.N., The Quantum-Field Renormalization Group in Critical Behavior Theory and 
Stochastic Dynamics [in Russian], PIYaF, St. Petersburg (1998); English transl.: The Field 
Theoretic Renormalization Group in Critical Behavior Theory and Stochastic Dynamics, 
Chapman and Hall/CRC, Boca Raton, Fla. (2004). 

[28] Karl F. Freed, Adv. Chem, Phys., 22, p. 1 (1972). 

[29] M.V.Komarova, M.Yu.Nalimov, Theoretical and Mathematical Physics, 126, p. 339 (2001). 



